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ABSTRACT 


A  formulation  of  th.2  Boltzmann  equation  for  dc  and  microwave  electric 
field  excitation  in  a  uniform,  multicomponent  gas  has  been  expressed  in  a 
computer  program  at  the  National  Bureau  of  Standards  by  L.  R.  Megill.  This 
program  was  intended  for  use  in  the  study  of  atmospheric  afterglow  by  NBS 
but  was  given  to  the  Aerospace  Corporation  for  conversion  tc  study  microwave 
breakdown.  The  conversion  of  the  program,  the  alteration  of  the  collision 
cross-section  data  pack,  and  the  numerical  verification  that  the  program  plus 
data  pack  accurately  describe  the  response  of  air  to  high  intensity  microwaves 
are  described  in  this  report. 
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I.  IHTROIXJCTOOK 


Prediction  of  the  electrical  breakdown  of  any  gas  requires  a  knowledge  of 
the  relationship  between  the  net  ionization  frequency  and  the  applied  electric 
field.  .ca  establishing  this  relationship  have  been  acquired  for  both  ac  and 
dc  fields  in  a  variety  of  gases.  Successful  theoretical  prediction  of  this 
relationship  has  been  restricted  to  the  rare  gases  because  of  the  complexity 
of  the  Boltzmann  equation,  which  is  used  to  describe  the  behavior  of  the  gas 
under  the  inflv^nce  of  an  electric  field.  The  comutational  problems  associated 
with  a  complex  gas  having  many  energy  levels  may  be  avoided  for  the  rare  gases 
having  only  a  few  excitation  processes  in  addition  to  the  elastic  and  ionization 
collision  processes.  For  a  diatonic  gas  such  as  air  at  low  temperatures,  the 
multitude  of  rotational,  vibrational,  and  electronic  excitation  processes  greatly 
complicates  the  problem. 

However,  a  few  years  ago  Carleton  and  Megill  (Ref.  l)  at  the  National  Bureau 
of  Standards  (NBS)  developed  a  computer  program  for  studying  dry  air  under  atmos¬ 
pheric  conditions.  Their  studies  were  limited  to  electronic  state  excitation 
processes  that  required  field  strengths  well  helow  breakdown  levels.  The  appli¬ 
cation  of  their  program  to  breakdown  conditions  was  later  performed  by  Lombardini 
(Ref.  2).  His  report  was  limited  in  that  a  detailed  comparison  of  the  theoreti¬ 
cally  predicted  and  experimentally  observed  net  ionization  frequencies  was  not 
presented.  Such  a  comparison  would  be  valuable  as  an  aid  to  interpretation  of 
current  methods  of  presenting  reiined  experimental  results. 

The  Carleton  and  Megill  program  was  a  sound  analysis,  but  the  output  of  the 
program  would  benefit  considerably  from  improvement  of  the  electron  collision 
cross  sections  that  were  input  to  the  program.  Since  the  time  of  the  original 
compilation  of  the  data  pack,  the  accuracy  of  the  experimental  cross  sections 
has  improved  considerably.  Also,  the  role  of  Importance  of  the  different  in¬ 
elastic  processes  shifts  as  the  average  energy  of  the  electrons  increases.  For 
these  reasons,  the  data  pack  was  completely  updated.  Some  cross  sections  were 
added,  especially  in  the  electronic  excitation  energy  region,  which  is  of  prime 
importance  in  the  prebreakdown  gas  energy  distribution.  In  all  cases,  the  most 
recently  determined  or  highest  confidence  values  were  used. 


Solutions  otf  the  Boltzmann  equation  were  obtained  for  various  mixtures  of 
Ct,  and  llg  over  a  range  of  total  densities  with  ac  and  dc  driving  fields  ranging 
from,  below  breakdown  to  substantially  above  breakdown.  The  principal  results 
axe:  the  electron  distribution  function,  average  energy,  mobility,  diffusion 
coefficient,  and  excitation  frequencies  for  the  various  molecular  states  including 
the  total  collision,  ionization,  and  attachment  frequencies. 
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n.  THE  BCUrZMftHI  EQUATION  solved 


The  derivation  of  the  Boltzmann  equation  is  covered  in  detail  in  the 
original  Carleton  and  Me  gill  (Ref.  1)  article  and  a  variety  of  other  sources 
(Ref.  2)  and  will  not  he  presented  here.  The  derivation  starts  with  the  Bolts- 
nann  equation  for  spatially  uniform  electron  distributions  t(y,  t) 


df 

dt 


+  a 


(V/) 


(i) 


coll 


(1) 


where  a  is_the  acceleration  due  to  the  driving  force  and  Vy  is  the  gradient  in 
velocity  (v)  space.  The  solution  is  assumed  to  have  the  form 

rft  s  f0(?)  *  e*  •  "v  f^v)  +  (b  x  e)  •  v  f2(v) 


where  E  and  B  are  the  electric  and  magnetic  fields,  respectively.  An  expansion 
of  f(v)  in  spherical  harmonics  in  "v  space  results  in  the  following  equation: 


V  A  fa  m  H 

V  (eE)2a(u)k  M1  V  E 


12kW  o’  (u) 

<Le(u)  g  (u)  +  - — -  g(u) 

in  sfa 


7®  E/  [<u  *  V  *  uu>  «<u  *  uij>  -  nj(u)e(u)]du}  <2' 


where  E  is  the  nns  amplitude  of  the  electric  field  with  angular  frequency  m, 
g(u)  -  fQ(u);  u(eV)  =  kv  ,  where  k  is  the  units  conversion  factor;  (e/m)  is  the 
electron  charge  to  mass  ratio;  is  the  density  of  species  i  with  mass  M^,  which 
has  an  elastic  cross  pection  tric(u)>  8X1  average  rotational  cross  section  of  <r  , 
and  inelastic  processes  j  described  by  °ij(u)  with  threshold  energies  u^eV) . 

is  the  average  energy  loss  in  a  rotational  collision  process  with  species  i. 
The  expressions  for  f  (u)  and  f  (u)  are  given  in  terms  of  frt(u)  or  g(u)  as 

1  O  v 
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f,<»)  •  -  2Z  a(«)  g 
’  $T  fl(u) 

where  Q  m  (eB/m), 

v(n)  (a,2  +  n2  +  */2(u)) 

e-W  *  —5 - 5 - r - 5 - 5 - 5- 

(r  -  if)  ♦  V  (u)  ♦  2i/(u)  («2  ♦  if) 


(3) 

(*) 

(5) 


(6) 


The  valid  solutions  for  large  field  strengths  axe  limited  to  small  fractions  of 
fj/f0  and  tg/f q. 

In  the  driving  frequency  domain,  either  the  solutions  are  limited  to 

_1 
U) 

V 


0>lmm 

“/  , 


du 


6(u)u*^u) 


=  relaxation  time 


(7) 


where  ®  (u)  is  the  average  fractional  energy  loss  of  an  electron  per  collision 
at  energy  u,  or  w  may  be  zero.  Thus,  there  is  a  low  driving  frequency  region 
that  would  produce  solutions  of  questionable  validity.  An  alternate  expression 
(Ref.  3)  requires  that 


e 

m 


A 

be  small.  Roughly,  this  expression  means  that  (e/m)(E/u,)  be  small, 
gradients  are  negligible  provided 


1  >  2E  —  + 
<u> 


X 

A 


Spatial 
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where  X  is  the  electron  mean  free  path,  <u>  is  the  overage  electron  energy,  and 
A  is  the  diffusion  length. 

Further  discussion  of  the  Units  of  the  theory  In  relation  to  observable 
quantities  may  be  found  in  McDonald's  (Ref.  3)  consideration  of  the  Units  of 
diffusion  theory. 


V 
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m.  CALIBRATION  QF  THE  COMPUTER  PROGRAM 


Frequently,  computer  programs  are  written  and  used  on  complex  problems 
without  any  attempt  to  solve  a  classical  problem  that  is  analytically  tractable. 
Such  an  exactly  solvable  problem  should  always  be  the  first  result  of  the  program 
In  order  to  establish  a  reference  of  accuracy  for  the  solution  mechanism. 

The  Maxwellian  distribution  is  an  exact  solution  of  the  Boltzmann  equation 
with  the  collision  frequency  equal  to  a  constant.  This  problem  was  solved  by 
the  routine  and  produced  numerical  representations  of  the  Maxwellian  accurate  to 
four  places  as  checked  by  use  of  a  standard  four-place  table  (Ref.  4);  the  actual 
solution  may  be  many  places  more  accurate. 

A  solution  similar  to  a  Maxwellian,  but  with  an  average  energy  given  by 
kT  +  (const)  E  where  k  is  the  Boltzmann  constant  and  T  ic  temperature,  results 
when  the  constant  collision  frequency  problem  has  a  dc  electric  field  of  strength 
E  applied.  This  solution  was  also  obtained  to  at  least  four  places. 

Further  calibration  of  the  program  to  include  the  data  pack  will  be  presented 
in  Section  VI. 
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IV.  IONISATION 


In  ‘breakdown  studies,  the  production  of  free  electrons  by  ionization  is 
the  single  most  important  electron  production  mechanism.  The  original  NBS  pro¬ 
gram  and  a  similar  program  developed  at  Westinghouse  by  Phelps  treat  the  ioniza¬ 
tion  process  the  same  as  any  inelastic  excitation  process  in  that  the  primary 
electron  is  considered  to  have  lost  an  amount  of  energy  equal  to  the  threshold 
energy  and  the  secondary  electron  is  ignored. 

Actually,  the  two  electrons  resulting  from  the  ionization  are  free,  and 
they  share  the  energy  of  the  primary  electron.  The  two  resultant  electrons  are 
indistinguishable  and  individually  may  have  energies  ranging  fran  zero  to  E  -  ^ 
(E  =*  energy  of  the  primary  electron).  In  the  Aerospace  version  of  the  program, 
it  is  assumed  that  the  resultant  electrons  share  the  available  energy.  That  is, 
they  each  reenter  the  electron  distribution  at  an  energy  §-(E  -  E^).  This  choice 
represents  an  average  in  the  available  energy  distribution  between  the  two  elec¬ 
trons  for  the  case  of  reasonably  low  impact  energies  as  encountered  in  the  break¬ 
down  environment.  Or,  it  may  be  stated  that  it  is  highly  unlikely  that  one 
electron  acquires  all  the  available  energy  in  the  ionization  process  and  the 
other  electron  none. 
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V.  CROSS -SECTION  DA TA 


The  cold-air  data  pack  cross  sections  in  Og  and  Ng  over  the  energy  range 
from  0  -  20  eV  are  those  described  In  Ref.  5  but  edited  to  eliminate  processes 
of  minimum  Importance  to  the  breakdown  problem.  The  cross  sections  used  In  the 
data  pack  are  listed  In  Table  I  In  their  order  of  occurrence  in  the  program. 

For  the  most  part,  the  determinations  of  Phelps  (Ref.  6),  Schulz  (Ref.  lt, 

Chen  (Ref.  8),  and  Rapp  and  Brlglla  (Ref.  9)  have  been  relied  upon  heavily.  The 
exception  is  the  Og  clastic  cross  section,  which  is  obtained  by  transformation 
properties  from  Ng  as  described  in  Ref.  5« 

The  best  confidence  that  can  be  ascribed  to  some  of  the  cross  sections  is 
±5$,  and  the  worst  is  ±20$.  These  confidences  are  the  estimates  of  Phelps 
(Ref.  6)  and  are,  to  a  certain  extent,  arbitrary.  For  the  study  of  breakdown 
and  prebreakdown,  the  most  arbitrary  portion  of  the  data  pack  is  comprised  of 
the  ’’lumped  electronic"  excitation  cross  sections.  This  is  due  to  the  fact 
that  the  electronic  states  that  are  "lumped"  lie  in  just  the  energy  range  forming 
the  last  "energy  loss  barrier"  or  energy  loss  mechanism,  which  keeps  the  tail 
of  the  electron  distribution  from  rising  in  value  in  the  energy  range  of  the 
ionization  process.  The  rising  in  the  tail  of  the  distribution  function  results 
in  an  increase  in  the  ionization  rate. 
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Table  I.  Air  Data  Pack  Cross-Section  Composition 


Species 

ID  No. 

Cross  Section 

Threshold 

(eV) 

Energy  Loss 
(u=energy) 

°2 

1 

Elastic 

0 

0 

2 

Elastic 

0 

N2 

3 

Elastic 

0 

a"eu'Y 

Black 

4 

d 

N2 

5 

Rotational 

2.49  x  10"4 

Blank 

6 

N2 

7 

X  -  c\ 

11.2 

11.2 

°2 

8 

vlb 

0.196 

0.196 

Blank 

9 

Blank 

10 

Blank 

11 

°2 

12 

X  -  a1* 

O.98 

0.98 

4 

13 

I.65 

1.65 

°2 

14 

x-a^4 

4.47 

4.47 

°2 

15 

X-B^r 

6.17 

6.17 

92 

l6 

Lumped,  elect 

9.5 

9-5 

°2 

17 

Dlssoc  att 

4.0 

4.0 

°2 

18 

Ionization 

12.2 

cvj 

1 

U"\ 

• 

O 

0 

19 

3p  - 

I.96 

1.96 

0 

20 

3P  -  Xs 

4.17 

4.17 

0 

21 

Lumped  elect 

9-5 

9.5 

0 

22 

Ionization 

13.6 

0.5(u  -  13.6) 

n2 

23 

l»  «  0  -  1 

0.29 

0.29 

2k 

0-2 

1.70 

1.70 

"a 

25 

O 

1 

UJ 

1.80 

1.80 

"2 

26 

0-4 

2.0 

2.0 

aFrcra  Ref.  5. 
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Table  I.  Air  Data  Pack  Cross-Section  Composition  (Concluded) 


Species 

ID  Ho. 

Cross  Section 

Threshold 

(eV) 

Energy  Loss 
(u=energy) 

N2 

27 

o 

I 

VJI 

2.0 

2.0 

n2 

28 

o  -  6 

2.2 

2.2 

"a 

29 

0-7 

2.4 

2.4 

N2 

30 

0-8 

2.6 

2.6 

h 

N2 

n2 

31 

32 

33 

X  - 

X-B3/ 

X  -  a1- 

g 

6.7 

7.39 

8.4 

6.7 

7.39 

8.4 

n2 

3* 

Lurped  elect 

5.0 

5.0 

n2 

35 

Phelps  14  eV 

14. 

14. 

n2 

36 

Ionization 

15-58 

0.5(u  -  15.58) 
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VI.  DC  CALIBRATION 


The  final  data  pack,  as  described  in  Table  I  and  in  Section  V,  was  used  in 
a  calibration  sequence  to  establish  confidence  in  future  calculations.  The 
standard  for  calibration  was  Phelps’  N2  results  in  a  dc  field.  Phelps’  (Ref.  6) 
results  closely  follow  the  experimental  drift  tube  measurements,  thus  giving  the 
present  calibration  a  reference  to  experimental  values. 

Because  the  underlying  interest  in  this  study  is  in  breakdown  and  prebreak¬ 
down  field  strengths,  the  calibration  to  Phelps’  (Ref.  6)  results  was  limited  to 

this  regime.  The  test  conditions  were  at  a  temperature  of  T7°K  with  E/N  values 

-IS  p 

of  1,  2,  3,  4,  and  5  x  10  7  V-cm  .  Comparison  was  made  with  Phelps’  (Ref.  6) 

curve  of  ionization  coefficient  a^/N  vs  e/n,  where  N  is  the  gas  density  (Ref.  6, 
Fig.  ?),  and  with  his  distribution  function  for  e/n  =  3  x  10-1^  V-cm2  (Ref.  6, 
Fig.  8). 

Calculation  of  *  t^w,  where  is  the  ionization  frequency  and  w  is  the 
drift  velocity,  was  made  for  a  range  of  e/n  values .  Values  of  Phelpa '  character- 


The  agreement  was  within  the  same  error  as  Phelps'  results  with  respect  to  the 
experimental  values  and  was  as  close  as  was  readable  from  the  graphs  (Figs.  1 
and  2).  Similarly,  the  distribution  function  for  the  case  e/n  =  3  X  10**^  V-cm2 
was  as  close  to  Phelps’  results  (Ref.  6,  Fig.  8)  as  could  be  read  from  the  graph. 

At  this  point,  it  must  be  noted  that  it  was  necessary  to  extend  the  range 
of  integration,  insofar  as  the  effects  of  ionization  and  elastic  processes  were 
concerned,  to  40  eV  for  the  higher  values  of  e/n.  This  was  necessary  to  obtain 
a  meaningful  approximation  to  the  ionization  rate  integral,  because  the  tail  of 
the  distribution  function  between  20  eV  and  30  to  35  eV  gave  a  nonnegligible 
contribution  to  the  rate  for  higher  values  of  e/n.  For  the  case  of  ac  fields, 
the  problem  of  the  high  energy  tail  of  the  distribution  function  does  not  develop 
until  fields  considerably  above  breakdown  are  encountered.  Also,  it  was  noted 
that  the  f^(u)  term  of  the  expansion  became  a  large  fraction  of  fQ  for  e/n  greater 
than  4  x  lcf1^  V-cm2. 
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E/N,  V-cm2 


Plg.l.  Comparison  of  drift  velocity  and  Phelps*  characteristic 
energy  for  Ng  in  a  dc  field  in  the  K/N  range  of  interest 
in  breakdown  studies.  [  x  *s  represent  the  results  of 
typical  calculations  performed  in  this  study;  other  sym¬ 
bols  represent  data  from  Ref.  6.  Illustration  is  from 
Phelps  (Ref .6,  Fig A),  courtesy  of  the  American  Insti¬ 
tute  of  Physics.] 


Fig. 2.  Comparison  of  the  first  Townsend  ionization  coefficient 
for  Np  in  a  dc  field  in  the  e/n  range  of  interest  in 
breakdown  studies.  [  x  's  represent  the  results  of 
typical  c  ilculations  performed  in  this  study;  other 
symbols  Tepresent  data  from  Ref.  6.  Illustration  is 
from  Phelps  (Ref.  6,  Fig. 7), courtesy  of  the  American 
Institute  of  Physics.] 


Figure  1  shows  the  reported  values  of  Phelps’  characteristic  energies 
and  the  drift  velocity  w  with  the  calibration  run  values  shown  as  X's.  Figure 
2  shows  the  agreement  of  the  Townsend  Ionization  coefficient;  the  computational 
results  are  shown  as  X's.  The  disagreement  in  this  figure  becomes  appreciable 
for  values  of  E/N  >  3  X  10  V-cm  .  At  this  point,  the  first  order  expansion 
term  Is  about  0.75  of  the  zeroth  order  term,  so  the  spherical  distribution  function 
approximation  Is  far  from  valid.  Also,  about  10$  of  the  electron  energy  loss  is 
through  ionization,  which  means  that  new  electrons  are  being  generated  at  a  non- 
negligfble  rate.  These  conditions  are  well  above  breakdown  and,  hence,  will  not 
effect  the  eventual  purpose  of  the  program. 

A  further  comparison  of  the  difference  between  Phelps'  solution  and  the 
present  one  Is  provided  at  e/n  »  3  X  10 V-cm2  where  he  shows  a  value  of 
is  4.6  eV  and  the  present  result  Is  4.3  eV.  This  difference  may  come  from  the 
ionization  process  treatment  discussed  In  Section  IV. 

Further  close  agreement  Is  found  with  the  average  energy  defined  as  1.5 
times  the  diffusion  coefficient  divided  by  the  mobility  as  shown  in  Brown  (Ref. 

10).  This  particular  plot  points  out  the  fallacy  in  the  assumption  by  MacDonald 
(Ref.  3>  P*  176)  that  the  average  electron  energy  is  linear  in  e/p,  where  p  is 
the  pressure,  when  the  experimental  data  show  (Ref.  10)  that  the  variation  is 
logarithmic  in  E/p. 
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VII.  HIGH-PRESSURE  MJCROHAVB  KIPHUMMEAL  COKRARISORf 


FOr  a  driving  frequency  of  9*375  GHz  In  air  at  100  torr,  the  electron 
collision  frequency  is  much  less  than  the  driving  frequency,  and  the  "effective 
field"  frequently  used  to  correlate  experimental  data  is  practically  the  same 
as  the  applied  field.  This  fact  was  utilized  to  provide  a  further  check  of  the 
program's  output  by  comparison  with  experimental  plots  of  net  ionization  frequency 
vs  "effective”  electric  field  strength. 

Thompson  has  verified  the  ionization  rate  data  obtained  by  other  researchers 
(Ref.  11).  Thompson's  (Ref.  12)  results  are  presented  in  Fig.  3  with  the  calcu¬ 
lated  results  obtained  from  the  computer  program.  The  agreement  is  seen  to  be 
qualitatively  and  quantitatively  correct  and  falls  veil  within  the  experimental 
error  bars. 

For  the  9-375  GHz  field  necessary  to  breakdown  air  at  100  torr,  the  value 
calculated  by  the  program  is  between  3000  and  3500  V«cm.  This  value  lies  on  the 
curve  given  by  *fec Donald  (Ref.  3,  p.  165). 

6  2 

The  mobility  calculated  at  breakdown  is  1.44  x  10  V-cm  -  sec,  which  may 

be  compared  with  a  result  from  an  analysis  attributed  to  Brown  (Ref.  3,  p.  172), 

6  2 

which  gives  a  value  of  1.25  X  10  V-cm  -eec.  This  is  a  disagreement  ox'  less  than 
15#  with  the  experimental  mobility. 

Agreement  with  low-pressure  "effective  field"  results  is  not  forthcoming 

because  the  "effective  field"  concept  is  not  valid  because  of  the  strong  variation 

of  collision  frequency  with  electric  lie  Id  strength  as  nay  be  seen  in  Fig.  4. 

It  is  not  too  surprising  to  see  chat  the  collision  frequency  exhibits  a  dependence 

on  electric  field  strength  and  frequency  as  well  as  pressure.  An  examination  of 

Eq.  (2)  for  the  time -dependent  case  shows  that  v,he  differential  equation  cannot 

transform  in  such  a  way  as  to  exhibit  the  "effective  field"  results  unless  certain 

conditions  are  met.  For  the  net  (1/.  -  v  )/p  to  be  invariant  with 

1  & 
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Fig. 3*  Net  ionization  frequency  for  air  calculated  with  an  ac 
driving  field  of  9»^kMc  in  comparison  with  experimental 
values  from  Refs.  11  and  12,  plotted  as  a  function  of 
effective  field  for  various  pressures.  [illustration 
is  from  Epstein  and  Lenander  (Ref.  13,  Fig.  3).] 


Fig.  4.  Average  total  collision  frequencies  as  a  function  of 
effective  field  for  various  pressures,,  shown  with  the 
frequently  suggested  value  of  5*3  x  10  .  [illustration 
is  from  Lenander  (Ref. 5).] 


the  "effective  field"  divided  by  pressure  oust  be  a  constant.  This  is  required 

so  that  the  electron  distribution  function  shape  will  be  Independent  of  pressure 

and  produce  the  same  collision  frequency  integral.  This,  in  turn,  is  obtainable, 

consistent  with  the  assumption  that  v  varies  linearly  with  pressure,  only  if 

c 

the  energy-dependent  total  collision  frequency  is  replaced  by  a  constant.  If  v 

V 

is  energy  dependent,  as  it  must  be,  the  nonlinearity  of  Eq.  (2)  emphasises  the 
error  inherent  in  the  "effective  Held"  approach. 

Free  these  caoaparlscns,  it  appears  that  the  calculated  observable  quantities 
for  9-37*  GHz  excitation  at  p  «  100  ton*  agree  quite  veil  with  experimental  values. 
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VUI.  SUMMATION  AND  CONCLUSIONS 


The  KBS  computer  program  for  solution  of  the  Boltzmann  equation  using  a 
multicross -section  model  for  air  was  converted  and  updated  for  use  at  Aerospace 
for  prediction  of  microwave  Breakdown  This  involved  a  reevaluation  of  the  Basic 
theory,  a  more  physical  treatment  of  the  ionization  products,  and  a  determination 
of  the  most  reliable  and  representative  cross  sections  in  air. 

The  program  was  tested  for  determination  of  its  ability  to  solve  the  Boltz¬ 
mann  equation  exactly  for  the  Maxwell  solution.  The  data  pack  of  cross  sections 
was  verified  as  Being  representative  of  air  By  calculation  of  the  observable 
quantities:  electron  mobility,  diffusion  coefficient,  ionization  coefficient, 
and  net  ionization  frequency.  Good  agreement  was  obtained  with  values  already 
known  from  experiment. 

On  the  Basis  of  the  calculated  test  results,  it  is  felt  that  the  program 
and  the  data  pack  will  produce  meaningful  predictions  of  microwave  Breakdown  in 
air. 
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APPENDIX 

The  Computer  Program 

The  program  produces  an  inward  integration  of  Eq.  (2)  utilizing  a  data 
table  of  cross  sections  defined  on  each  0.10  eV  from  0  to  20.0.  The  data 
table  is  in  the  form  of  a  binary  tape,  which  has  the  cross-section  thresholds 
as  the  first  record  and  the  36  cross  sections  (some  are  zero  dumnys)  as  the 
next  36  records. 

Data  input  directly  to  the  program  consists  of  the  density  fractions  of 
°2>  °>  and  amu  masses  of  Og,  0,  and  Ng  (must  not  be  zero);  gas  temperature,  °K; 
pressure,  torr;  rms  electric  field  strength,  V-cm;  driving  frequency,  cps; 
magnetic  field  perpendicular  to  the  electric  field,  gauss;  and  the  maximum 
energy  for  the  integration,  eV(£20eV).  The  units  of  the  integration  are  a  mix¬ 
ture  of  esu,  cgs,  and  eV. 

The  energy  step  of  the  integration  is  the  energy  range  divided  by  2000. 

As  the  inward  integration  proceeds,  the  necessary  cross-section  values  are  ob¬ 
tained  by  linear  interpolation  of  the  data  table  values  at  the  nearest  0.1  eV 
points . 

As  the  inward  integration  proceeds,  the  distribution  function  value  is 
tested  to  prevent  overflow,  which  can  result  from  too  large  an  energy  range, 
and  the  function  derivative  is  tested  to  find  the  energy  at  which  the  Lorentz 
term  is  the  largest.  Step  by  step,  the  observables  or  moments  of  the  distribution 
function  are  calculated.  The  proper  normalization  is  supplied  at  the  end  of  the 
integration. 

The  output  quantities  are:  the  gas  constituent  densities;  gas  temperature; 
total  gas  pressure;  energy  increment  for  the  integration  step;  the  driving  electric 
field  and  its  frequency;  the  magnetic  field  strength;  the  effective  field  based  on 
the  calculated  total  collision  frequency;  the  average  electron  energy;  average 
electron  speed;  total  collision  frequency;  conductivities  parallel  and  perpendi¬ 
cular  to  the  electric  field  and  the  Hall  magnetic  conductivities  per  electron; 
diffusion  coefficient;  drift  velocity;  Phelps'  characteristic  energy  defined  as 
the  diffusion  coefficient  divided  by  the  drift  velocity,  expressed  in  eV;  mean 
free  path;  average  interparticle  distance;  a  table  of  the  rate  constants  for 
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ea-'h  ot  the  cross  sections  considered  and  the  amount  of  electron  energy  loss 
to  each;  Q,  radiative  attachment  rate;  0^  dissociative  attachment  rate;  total 
ionization  rate;  the  normalization  factor;  the  largest  ratio  of  the  Larentz  term 
to  the  zeroth  order  term  and  the  energy  at  which  it  occurs;  and  a  listing  of 
this  distribution  function  at  each  tenth  point  of  the  integration  range  on  the 
corresponding  energy.  All  the  above  quantities  are  written  out  in  cgs-esu  or 
are  clearly  labeled. 

When  the  calculations  for  one  case  are  complete,  the  program  looks  for  a 
new  set  of  input  data  before  terminating.  No  provisions  are  made  within  the 
program  that  prevent  improper  solutions  or  indefinite  numbers.  It  is  assumed 
that  the  data  input  to  the  program  is  correct  and  of  a  reasonable  nature  such 
that  problems  of  that  type  will  not  arise. 
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